THE MOMENTS OF MINKOWSKI QUESTION MARK FUNCTION: THE DYADIC 

PERIOD FUNCTION 



GIEDRIUS ALKAUSKAS 

Abstract. The Minkowski question mark function l(x) arises as a real distribution of rationals in the 
Farey tree. We examine the generating function of moments of ?(x). It appears that the generating 
function is a direct dyadic analogue of period functions for Maass wave forms and it is defined in the cut 
plane C \ (0, oo). The exponential generating function satisfies the integral equation with kernel being the 
Bessel function. The solution of this integral equation leads to the definition of dyadic eigenfunctions, 
arising from a certain Hilbert-Schmidt operator. Finally, we describe /5-adic distribution of rationals in 
the Stern-Brocot tree. Surprisingly, the Eisenstein series G\(z) does manifest in both real and /?-adic 
cases. 
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1. Introduction 

This papeiQ is the first in the series of four papers (others being JTJ, [0 and (3]|) which are devoted 
to the study of moments and integral transforms of the Minkowski question mark function. 

The function ?(jc) ("the question mark function") was introduced by Minkowski in 1904 [|34l as 
an example of continuous and monotone function ? : [0, 1] — > [0, 1], which maps rationals to dyadic 
rationals, and quadratic irrationals to non-dyadic rationals. It is though more convenient to work 
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Figure 1. The Minkowski question mark function F(x), x e [0, 2] 

with the function F(x) :=?(-^f), x e [0, oo) u {oo}. Thus, for non-negative real x it is defined by the 
expression 



F([a ,a u a 2 ,a 3 ,...]) = 1-2^+2 



(ao+ai) _ 2"(«o+ai+a2) 



+ 



(1) 



where x = [ciq, Cl\, Cl2i ■■■ 

] stands for the representation of x by a (regular) continued fraction ll23Tl . 
Hence, according to our convention, ?(x) = 2F(x) for x e [0, 1]. For rational x the series terminates 
at the last nonzero element a n of the continued fraction. 

1.1. Short literature overview. The Minkowski question mark function was investigated by many 
authors. In this subsection we give an overview of available literature. 

Denjoy [fTTI gave an explicit expression for F{x) in terms of continued fraction expansion; that is, 
formula CD. He also showed that ?(x) is purely singular: the derivative, in terms of the Lebesgue 
measure, vanishes almost everywhere. Salem [|43l proved that ?(x) satisfies the Lipschitz condition of 
order 5^-, where y = ^-y^, and this is in fact the best possible exponent for the Lipschitz condition. 

The Fourier-Stieltjes coefficients of ?(x), defined as J* e 2mnx d?(x), where also investigated in the 
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same paper (these coefficients also appeared in [6]; see also IflTTO . The author, as an application 
of Wiener's theorem about Fourier series, gives average results on these coefficients without giving 
an answer to yet unsolved problem whether these coefficients vanish, as n — » oo (it is worth noting 
that [1] analogous Fourier coefficients are introduced and examined). Kinney [|24l proved that the 
Hausdorff dimension of growth points of 7(x) (denote this set by Jl) is equal to a = |( J Q log 2 (l + 
x) dl(x)) 1 . Numerical estimates for this constant were obtained in E71 and 11481 : based on the three 
term functional equation, we are able to calculate the Kinney's constant to a high precision (in the 
appendix of [3] we calculate 35 exact digits; note that some digits presented in |[38ll are incorrect). 
Also, if xq € ?(jc) at a point xq satisfies the Lipschitz condition with exponent a. The function ?(jc) 
is mentioned in IfTOll in connection with a "box" function. In [|28l Lagarias and Tresser introduce the 
so called Q-tree: an extension of the Farey tree which contains all (positive and negative) rationals. 
Tichy and Uitz [48] extended Kinney's approach (mainly, calculation of a Hausdorff dimension) to 
a parametrized class of singular functions related to ?(x). Bower considers the solution of the 
equation ?(x) = x, different from x = 0, \ or 1 . There are two of them (symmetric with respect to 
x = i), the first one is given by x = 0.42037233+ Q3J. Apparently, no closed form formula exists 
for it. In lfT2l Dilcher and Stolarsky introduced what they call Stern polynomials. The construction 
is analogous to similar constructions given in |fl"4l and Q. Nevertheless, in lfT2l all polynomials 
have coefficients and 1, and their structure is compatible with regular continued fraction algorithm. 
In [fT3l Dushistova and Moshchevitin find conditions in order T(x) = and = 00 to hold (for 
certain fixed positive real x) in terms of lim sup,^ ao+ai +- +a > an( j \[ m rmV^ " 0+a ' < K " + "' respectively, 
where x = \ao,a\,a2, ■■■] is represented by a continued fraction. The nature of singularity of ?(x) 
was clarified by Viader, Paradfs and Bibiloni 11381 . In particular, the existence of the derivative ?'(x) 
in R for fixed x forces it to vanish. Some other properties of ?(x) are demonstrated in lT3~9ll . In [|22l 
Kessebohmer and Stratmann studied various fractal geometric aspects of the Minkowski question 
mark function F(x). They showed that the unit interval can be written as the union of three sets: 
A := {x : F'(x) = 0}, A M := {x : F'(x) = 00}, and A^ := {x : F'{x) does not exist and F'(x) ± 00}. 
Their main result is that the Hausdorff dimensions of these sets are related in the following way: 
dim//(vf) < dim^(A^) = dim//(Aoo) = dim# («5f (h top )) < dim//(A ) = 1. Here Jf (h top ) refers to 
the level set of the Stern-Brocot multifractal decomposition at the topological entropy h top = log 2 of 
the Farey map Q, and dim#(vf- ) denotes the Hausdorff dimension of the measure of maximal entropy 
of the dynamical system associated with Q. The notions and technique were developed earlier by 
authors in ETI . The paper [26J deals with the interrelations among the additive continued fraction 
algorithm, the Farey tree, the Farey shift and the Minkowski question mark function. The motivation 
for the work 11371 is a fact that the function ?(jc) can be characterized as the unique homeomorphism 
of the real unit interval that conjugates the Farey map with the tent map. In [37J Panti constructs an 
n-dimensional analogue of the Minkowski function as the only homeomorphism of an n-simplex that 
conjugates the piecewise-fractional map associated to the Monkemeyer continued fraction algorithm 
with an appropriate tent map. In |[6l Bonanno and Isola introduce a class of 1 -dimensional maps 
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which can be used to generate the binary trees in different ways, and study their ergodic properties. 
This leads to studying some random processes (Markov chains and martingales) arising in a natural 
way in this context. In the course of the paper the authors also introduce a function p(x) =?(-fj), 
which is, of course, exactly F(x). Okamoto and Wunsch 113611 construct yet another generalization of 
?(jc), though their main concern is to introduce a new family of purely singular functions. Meanwhile, 
the paper by Grabner, Kirschenhofer and Tichy lfT8ll . out of all papers in the bibliography list, is 
the closest in spirit to the current article. In order to derive precise error bounds for the so called 
Garcia entropy of a certain measure, the authors consider the moments of the continuous and singular 
function F 2 ([cii, a 2 , ...]) = Z^i(-l)"~ 1 3~ (ai+ - +a,, ~ 1) (g„ + q n -\), where q+ stand for a corresponding 
denominator of the convergent to \a\, a 2 , ...]. Lamberger [29] showed that F(x) and F 2 (x) are the first 
two members of a family (indexed by natural numbers) of mutually singular measures, derived from 
the subtractive Euclidean algorithm. The latter two papers are very interesting and promising, and the 
author of this article does intend to generalize the results about F(x) to the whole family Fj(x), j e N. 

1.2. Stern-Brocot tree. Recently, Calkin and Wilf [8] (re-)defined a binary tree which is generated 
by the iteration 

a a a + b 

- i-> , , 

b a + b b 

starting from the root \ (this tree is a permutation of the well-known Stern-Brocot tree). Elementary 
considerations show that this tree contains any positive rational number once and only once, each 
represented in lowest terms [H. First four iterations lead to 

1 

1 



3 
1 

\ 

4 
1 

Thus, the nth generation consists of 2"" 1 positive rationals x„ , 1 < i < 2 n ~ l . We denote this tree by 7~, 
and its nth generation by T {n) . The limitation of this tree to the interval [0, 1] is the well-known Farey 
tree (albeit with different order). Reading the tree line by line, this enumeration of Q + starts with 
112132314352534 
T'2'T'3'2'3'T'4'3'5'2'5'3'4'T'""" 
This sequence was already investigated by Stern 11461 . where one encounters the definition of the 
Stern-Brocot tree. The sequence satisfies the remarkable iteration discovered by M. Newman [|35l : 

*i = 1> Xn+l = l/(2[x„] + 1 - x n ), 

thus giving an example of a simple recurrence which produces all positive rationals, and answering 
affirmatively to a question by D.E. Knuth. The nth generation of T consists of exactly those rational 




MOMENTS OF THE QUESTION MARK FUNCTION 



5 



numbers whose elements of the continued fraction sum up to n; this observation is due to Stern 
Il46ll . Indeed, this can be easily inherited directly from the definition. First, if rational number f is 
represented as a continued fraction [ao,«i, ...,a r L then the map | — » ^ maps | to [ao + \,a\..., a r ]. 
Second, the map | — > -jr maps | to [0, a { + l,...,a r ] in case | < 1, and to [l,a ,ai, ...,a r ] in case 
| > 1 . This simple fact is of utmost importance in our work: though it is not used in explicit form, this 
highly motivates the investigations of moments M L and m £ , given by ©. The sequence of numerators 
of the Calkin-Wilf tree 

0,1,1,2,1,3,2,3,1,4,3,5,2,5,3,4,1,... 

is called the Stern diatomic sequence [|46ll . OTTl and it satisfies the recurrence relations 

s(0) = 0, = 1, s(2n) = s(n), s(2n + 1) = s(n) + s(n + 1). (2) 

In the next section we will show that each generation of the Calkin-Wilf tree possesses a distribution 
function F n (x), and F n (x) converges uniformly to F(x). This is by far not a new fact. Nevertheless, 
we include the short proof of it for the sake of self-containedness. The function F(x), as a distribution 
function, is uniquely determined by the functional equation ©. This implies the explicit expression 
dU) and the symmetry property F(x) + F(l/x) = 1 . The mean value of F(x) was investigated by several 
authors ( ll4"7Tl . lt5T10 and was proved to be 3/2. We will obtain this result using quite a different method. 

1 .3. Motivation and description of results. The aim of this paper is to give a different treatment of 
Minkowski's All papers so far were concerned with F(x) as the function per se. Nevertheless, it 
appears that there exist several natural integral transforms of F(x), which are analytic functions and 
which encode certain substantial (in fact, all) information about the question mark function. Each of 
these transforms is characterized by regularity properties and a functional equation. Lastly and most 
importantly, let us point out that, surprisingly, there are striking similarities and analogies between 
the results proved here as well as in [Q]|, with Lewis' -Zagier's 11331 results on period functions for 
Maass wave forms. That work is an expanded and clarified exposition of an earlier paper by Lewis 
11321 . The concise exposition of these objects, their properties and relations to Selberg zeta function 
can be found in 11521 . The reader is strongly urged to compare results in this work with those in OBI . 
Let, for example, u(z) be a Maass wave form for PSL 2 (Z) with spectral parameter s. The similarity 
arises due to the fact that the limit value of u(z) on the real line, given by u(x + iy) ~ y l ~ s U(x) +y s U(x) 
as y — » 0+, satisfies (formal) functional equations U(x + 1) = U(x) and \x\ 2s ~ 2 U( - -) = U(x). Thus, 
these are completely analogous to the functional equations for F(x), save the fact that U(x) is only a 
formal function - it is a distribution (e.g. a continuous functional in properly defined space of func- 
tions). Thus, our objects G(z), m(f) and M{t) are analogues of objects i//(z), g(w) and (p(w) respectively 
(see Section 2 of this work and [33]). In HI it is shown that in fact L-functions attached to Maass 
wave forms also do have an analogue in context of the Minkowski question mark function. 

This work is organized as follows. In Section 2 we demonstrate some elementary properties of 
the distribution function F(x). Since the existence of all moments is guaranteed by the exponential 



6 



GIEDRIUS ALKAUSKAS 



decay of the tail, our main object is the generating function of moments, denoted by G(z). In Section 
3 we prove two functional equations for G(z). In Section 4 we demonstrate the uniqueness of solu- 
tion of this functional equation (subject to regularity conditions). Surprisingly, the Eisenstein series 
Gi(z) appears on the stage. In Section 5 we prove the integral equation for the exponential generating 
function. In Section 6 a new class of functions emerging from eigenfunctions of a Hilbert-Schmidt op- 
erator (we call them "dyadic eigenfunctions") is introduced. These are dyadic analogues of functions 
discovered by Wirsing ll50ll in connection with the Gauss-Kuzmin-Levy problem. In Section 7 we 
describe p-adic distribution of rationals in the Calkin-Wilf tree. In the final section some concluding 
remarks are presented. 

2. Some properties of the distribution 

The following Proposition was proved by many authors in various forms, concerning (very related) 
Stern-Brocot, Farey or Calkin-Wilf trees, and this seems to be a well-known fact about a distribution 
of rationals in these trees. For the sake of completeness we present a short proof, since the functional 
equations for G(z) and m(t) (see Sections 3 and 5) heavily depend on the functional equation for F(x) 
and are in fact reformulations of these in different terms. 

Proposition 1. Let F n (x) denote the distribution function of the nth generation, i.e., 

F n (x) = 2 1 -"#{j : xf < x). 

Then uniformly F n (x) — > F(x). Thus, F(0) = 0, F(oo) = 1. Moreover, F(x) is continuous, monotone 
and singular, i.e., F'(x) = almost everywhere. 

Proof. Let x > 1. One half of the fractions in the n + 1-st generation do not exceed 1, and hence also 
do not exceed x. Further, 

a + b a 

< x <=> - < x- 1. 



Hence, 



Now assume < x < 1 . Then 



b ~ b 
2F n+1 (x) = F n (x-l) + l, n>\. 



a x 
< x <=> - < 



a + b b 1 - x 



Therefore, 



2F n+1 (x) = F n (-^—). 

x 1 - x' 

The distribution function F, defined in the formulation of the Theorem, satisfies the functional equa- 
tion 

2FU) . ( F %~ » + 1 ,3, 
I F(^) if 0<x<l. 
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For instance, the second identity is equivalent to 2F{-^) = F(t) for all positive t. If t = [b ,bi, ...], 
then j^j- = [0, l,b ,bi, ..] for t > 1, and = [0,bi + l,b 2 , ...] for t < 1, and the statement follows 
immediately. 

Now define 5 n (x) = F(X) In order to prove the uniform convergence F n — > F, it is sufficient 

to show that 

sup|<5„(*)l <T n . (4) 

;t>0 

It is easy to see that the assertion is true for n = 1. Now suppose the estimate is true for n. In view of 
the functional equation for both F n (x) and F{x), we have 

28 n+l (x) = 5 n {-^—) 
v l - x' 

for < x < 1, which gives sup < v<1 |<5„ +1 (x)| < 2~ n ~ x . Moreover, we have 

2£„ + i0c) = 6 n (x- 1) 

for x > 1, which yields the same bound for in the range x > 1. This proves ©. As it was noted, 
the singularity of F(x) was proved in IfTTTl and it follows from Khinchin results on metric properties 
of continued fractions. ■ 

Since F(x) has a tail of exponential decay (1 - F(x) = 0(2~ x ), as it is clear from (QQ)), all moments 
do exist. Let 

co oo 11 

M L = j x L dF(x), m L = J (-^-jfdFix) = 2 J x L dF(x) = J ^d?^). 



Therefore, M L and m L can also be defined as 

M L = lim2 1 -" V [a ,a u ..,a s ] L , m L = \im2 2 ~ n V [0, a u .., a s ] L , (5) 

n— *oo i i n — >co i i 

ao+ai+...+a s =n ai+,..+a s =n 

where the summation takes place over all rationals, whose elements of the continued fraction sum 
up to n. These expressions highly motivate our investigation of moments. Though the authors in 
[fT8ll considered the moments of F 2 (x) (see the introduction), it is surprising that the moments of 
Minkowski question mark function itself were never investigated. Numerically, one has 

Mi = 1.5, M 2 = 4.290926, M 3 = 18.556, M 4 = 107.03; 

mi = 0.5, m 2 = 0.290926, m 3 = 0.186389, m 4 = 0.126992. 

We will see that the generating function of m L possesses certain fascinating properties. Let u>(x) be a 
continuous function of at most polynomial growth: oj(x) = 0(x T ), x — » oo. The functional equation 
© gives F(x + n) = 1 - 2~ n + 2~ n F{x), x>0. Hence 

oo oo 1 1 oo 

f W (x)dF(x) = £ f co(x + n)dF(x + n)= f £ dF(x). 

n=0 " =0 
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Since, as noted above, F(x) has a tail of exponential decay, this integral does exist. Let x = t > 0. 
Since F(^) = \F(t), this change of variables gives 

^(777 + n) 



J (o(x)dF(x) = Y J J df(r) 



M=0 



(All changes of order of summation and integration are easily justifiable minding the condition on 

n 
2" 



o»(x)). Let o»(x) = x L , L e No. Then, if we denote 5^ = 2^=0 5J?r> we nave 



j>dF(*) = j^^^B^dFix). 



<=° 



Whence the relation 

L 



XmiV, ,L>0. (6) 
i=o ' 1 ' 

The exponential generating function of B L is 

5 ^ ~ XI 7T ? ~ 2j 2j 2^+>77 ~ ^ 2 n+1 ~ 2-e r 

L=0 ' L=0 n=0 ' n=0 

Denote by M(t) and m(t) the corresponding exponential generating functions of the coefficients M L 
and m L respectively. Accordingly, 

QO OO 1 

M(t) = J e xt dF(x), m(0 = J~ exp(-^-)dF(x) = 2 J" <?*dF(x). 


The relation © in terms of M(t) and m(£) reads as 

00 -MM -co - 

^ L! 2 - e ? f-t LI 2 - e' w V 

We see that the function m(f) is entire and M(t) has a positive radius of convergence. The last identity 
implies the asymptotic formula for M L . 



Proposition 2. For L e N , 



Proof. By Cauchy's formula, for any sufficiently small r, 

u r M( Z ) 
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Changing the path of integration, we get by the calculus of residues 

i(z) dz 



I m(z) \ L\ C m(; 



where R satisfies log 2 < R < | log 2 + 2ni\ (which means that there is exactly one simple pole of the 
integrand located in the interior of the circle \z\ = R). It is easily seen that the residue coincides with 
the main term in the formula of the Lemma; the error term follows from estimating the integral. ■ 

Also, © gives the inverse to linear equations ©: 

m L = M L -Y u M s \ , L>0. (8) 

Since B{t){2 - e') = 1, the coefficients B L can be calculated recursively: B L = J^Io (j^- Thus, 
B = 1, B } = 1, B 2 = 3,B 3 = 13, B 4 = 75, B 5 = 541. This sequence has number A000670 in Ii45l 
and traces its history back from Cayley. 

In the future we will consider integrals which involve m(t), and hence we need the evaluation of 
this function for negative t. 



Lemma 1. Let C = e' v lo § 2 = 0.4349+. Then C 2V? « m(-t) « C^' as t oo. 

exp(--^j-)dF(x). Hence, m(t) is positive for t e R. Let < M < 1. 
Since 1 - F{x) x 2~ x as x — » oo, and F{x) + F{\ /x) = 1, 

M oo 

m(-t) = ( J + J )exp(--^ T )dF(x)«2-^ + exp(- ] i^ T ). 

M 

This is valid for every M < 1 and a universal constant. A choice M = ^^ g2 gives the desired upper 
bound. To obtain the lower bound, note that 



M 

m(-t) > f exp(-^-) dF(x) » 2~ 1/M • exp ( - -^-). 
J x + 1 x M + 1' 



The same choice for M establishes the lower bound. Naturally, similar evaluation holds for the deriv- 
ative, since m'(-t) = X°° ITT ex P( _ lfr) dF(x). ■ 

We will prove one property of the function m(f) which represents the symmetry of F given by 
F{x) +F(l/x) = 1. 



Proposition 3. We have: m(t) = e'm(-t). 
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Proof. In fact, 

CO CO 

m(0 = j exp(-^ T )dF(x) = - J eX p(-^-)dF(l/x) 



CO CO 

= ( exp(— — )dF(x) = e' I exp( - -^—)dF(x) = m(-t)e' . 



o 

Whence the relations 



m L = W L )(-l) s m,, ,L>0. 



Thus, mi = m - m l5 which gives ni\ = 1/2, and this implies Mi = 3/2. Also, 2m 3 = -1/2 + 3m 2 . 
These linear relations are further investigated in [Q]|. 

3. The dyadic period function G(z) 
We introduce the generating power function of moments 

oo 

M(z) = Yj m lz l - 

L=0 

A priori, this series converges in the unit circle. Recall that J Q x n e~ x dx = T(n + 1) = n\. Thus, for 
real z < 1, the symmetry relation for m(f) gives: 

oo oo 

M( z ) = m(z0e"'dz 1 = J m(-z0e" ?(1 " z) d? = 
o o 

CO 

o 

Both integrals converge for z < 1 (since m L < 1, \m(z)\ < e z ), hence for these values of z we have the 
above identity. The function M(z) was initially defined for \z\ < 1; nevertheless, the above identity 
gives us holomorphic continuation of M(z) to the half plane %z < 1 /2. 



Lemma 2. The function Ai(z) can be analytically continued to the domain 



KX>\- 



J, CO 
exp(^t)dF(x). As noted above, |m(?)| < e for positive t (actually, 

Lemma \T\ combined with Proposition [3] gives a slightly better estimate). Therefore, for real z < 1 we 

have: 

CO oo CO 

M(z) = j j exp(-^zt)e-'dF(x)dt = j - ^ z dF(x). 
oo o x+1 

We already obtained the analytic continuation of M to the region {\z\ < 1} Ul^-z < 1 /2}. Let z = cr+jj 
with v > and o~ > 1/2. In the small neighborhood of z the imaginary part is bounded: y > yo > 0, 
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and also the real part is bounded: <x < cr . In this neighborhood the integral converges uniformly, 
because we have the estimate 



(1- 



x + 1 



z) 



> max) 



1 - 



x 



x+1 



-o~ 



x+ 1 



For < x < l/cr this gives the bound , and for x > 1/cr this gives the bound . Hence, in 
this neighborhood the function under the integral is uniformly bounded, which proves the uniform 
convergence of the integral and the statement of Lemma. ■ 



The system © gives us the expression of M L in terms of m s . In fact, there exists one more system 
which is independent of the distribution F(x); it simply encodes the relation among functions {j^) L 
and x s , given by 

A change x = gives 

And ultimately, 

^=z(lIiK (io) 

.v>L ' ' 

For the convenience, we introduce a function 

tx> 1 

G(z) ^Zl -_ -_ J^dFW =2/^^). (ID 

L_1 

Next theorem is our main result about G(z). The power series converges in the disc \z\ < 1 (including 
the boundary, as can be inherited from (flOl) ; moreover, this implies that there exist all left derivatives 
of G(z) at z = 1). The integral converges in the cut plane C \ (1, 00). 

Theorem 1. Let m L = j Q (-^) L dF(x). Then the generating power function, defined as G(z) = 
Zl=i m LZ L ~ l , has an analytic continuation to the domain C\R X>1 . It satisfies the functional equation 

~ t~~ 7T-^ G {^T—) + 2G ^ + V = W 
l-z (1 -zr v l -z' 

and also the symmetry property 

1/1 V 1 

G(z+l) = --G(- + l)--. 

z z y z ' z 

Moreover, G{z) = o(l) as z —* 00 and the distance from z to R.+ tends to infinity. 
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Proof. In analogy to M(z), for real z < define the following function: M Q (z) = J M(zt)e dt. 
In view of ©, this integral converges for real z < 0. Thus, 

CO CO CO 

M (z) = j j exp(xzt)e-'dF(x)dt= j j-^-dF(x). 



From argument akin to the one used in proving Lemma[2]we deduce that Mq(z) extends as an analytic 
function to the region C\R >0 . In this domain we see that 

M (z)-1 M(z+1)-1 



z z + 1 

which is the consequence of an algebraic identity 



(13) 



( 1 ) 1 - ( 1 

' 1 -XZ ' ~Z ~ ' 1 -r$7(z+ 1) ~ ' 



xz z 1_ 7TT^ +1 ) z+1 

The relation (fT3l) is independent of the specific distribution function, it simply encodes the information 
contained in (flOl) about the relation of x L to (x/(x + 1))'. On the other hand, the specific information 
about F(x) is encoded in © or (|7]). The comparison of these two relations gives the desired functional 
equation for G(z). In fact, for real t < the following estimate follows from © and Lemma CD 
\M(t)\ = |m(0(2 - eO _1 l ^ ^ 1; and thus for real z < we have: 



oo CO 

M( z ) = J m(zt)e-'dt= J (2 - e zt )M(zt)e~' dt = 







CO 

2M (z) - f M(z0e"' (1 " z) df = 2M (z) - ^0(7^)7^- 
J x l -z'l -z 



Finally, the substitution (1131) gives us the functional equation 

iz£ _ ^yw(-L) + 2^-M(z + 1) = M( Z ). 
1 +z 1 -z v l-z 7 z+1 

The principle of analytic continuation implies that this equation should be satisfied for all values 
of arguments in the region of holomorphicity of M(z). Direct inspection shows that for G(z) = )-1 
this equation reads as (fT2)) . Also, the symmetry property is a reformulation or ©. This proves the 
first part of the Theorem. 

Obviously, the last assertion follows from the integral representation of G(z) given by (TlTT) . ■ 

We call G(z) the dyadic period function, since its functional equation is completely parallel to a 
three term functional equations which are satisfied by rational period functions and period functions 
associated with Maass wave forms lT33ll . The word "dyadic" refers to the binary origin of the distribu- 
tion function F(x). Indeed, thorough inspection shows that the multiplier 2 in equations (1201) and (fT2~l) 
emerges exactly from the fact that every generation of T has twice as many members as a previous 
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generation. 

4. Uniqueness of G(z) 

In this section we prove the uniqueness of a function which have properties described in Theorem 
CD Note that two functional equations for G(z) can be merged into a single one. It is easy to check that 

- + \g(-) + 2G(z+1) = G(z) (14) 
z z z K z' 

is equivalent to both together. In fact, the change z h 1 /z in the equation (fl"4"l) gives the symmetry 
property, and application of it to the term G(l /z) of the above recovers the functional equation (fT2l) . 

Proposition 4. The function which satisfies the conditions ofTheoremUjis unique. 

Proof. Suppose there exist two such functions. Then their difference Gq(z) has the same behavior at 
infinity, and it satisfies the homogenic form of the equation (fT4l) . Let M = sup_ 1<t<0 |Gq(x)| = |Gq(xo)|, 
xq e [-1,0]. We will show that M = 0; by the principle of analytic continuation this will imply that 
Gq(z) = 0. Let z be real, -1 < z < 0. Let us substitute z i-» z - n in the equation (IT4b . n eM, n > 1, 
and divide it by 2". Thus, we obtain: 

Gp(z - n) _ G (z-n + 1) _ 1 / 1 \ n 

2" 2"" 1 2"(z-n) 2 %-n'" K } 

Note that for z in the interval [-1,0], ^ belongs to the same interval as well. Now sum this over 
n > 1 . The series on both sides are absolutely convergent, minding the behavior of G (z) at infinity. 
Therefore, 

°° 1 1 
-G Q (z) = V — -G ( ). 

2 n (z - nY y z-n> 

The evaluation of the right hand side gives: 

|Gofe)l * Z 2^ M = (l2 " 2 l0g22 ) M f ° r " 1 - Z - °" 

n=l 

The constant is < 1. Thus, unless M = 0, this is contradictory for z = x . This proves the Proposition. 
■ 

Note the similarity between (fl"5T) and the expression for the Gauss-Kuzmin-Wirsing operator W. 
The latter is defined for bounded smooth functions / : [0, 1] —* R by the formula 

oo 

The eigenvalue 1 corresponds to the function (see ll23l chapter III, for Kuzmin's treatment). The 
second largest eigenvalue -0.303663... (the Wirsing constant) leads to a function with no analytic 
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expression known [|50ll ; this eigenvalue determines the speed of convergence of iterates [W^/](x) 
to (for certain c e R). The spectral analysis of our operator is presented in Section 6. Paper [OQ 
contains much more details and results in this direction. 



Let 3z > 0. We remind that the Eisenstein series of weight 2 for PSL 2 (Z) is defined as p4| 

td ( m + nz ? 

(mind the order of summation, since the series is not absolutely convergent). This series has the 
following Fourier expansion: if q = e 2mz , then 

.2 



G l (z) = ^--Sn 2 Y j cr l (n)q". 



n=i 



Then this function is not completely modular, but we have the following identities ( Il44ll . chapter VII): 

G y {z + 1) = Gi(z), Gi(-l/z) = z 2 G x (z) - 2mz. 

Note that for 3z > 0, all arguments in (PT21) simultaneously belong to the upper half plane. It is 
surprising (but not coincidental) that the function j-G\(z) satisfies the functional equation (|T2|) for 
Sz > (see the remarks in Section 8 about possible connections in idelic setting). To check this 
statement, note that 

^ Gi ( ~ jry) = ~ 1)2Gi(z ~ 1} ~ 2ni(z ~ 1} ) = ~ z)2Gi(z) ~ (1 ~ z) - 



Thus, plugging this into (PT21) . we obtain an identity. If we define G\{z) = Gi(z) for 3z < 0, one 
checks directly that the symmetry property is also satisfied. This is a surprising phenomena. See the 
last section of 0]| for more speculations on this topic, where the space of dyadic period functions in 
the upper half plane (denoted by DPF°) is introduced. 

We end this section with presenting a system of linear equations which the moments m L do sat- 
isfy. This system is derived from the three therm functional equation (PT4l) and is a superior result in 
numerical calculations: whereas directly from the definition we can recover only a few digits of the 
moments, this method allows to calculate up to 60 digits and more. 

Proposition 5. Denote cl = Z^ti 2^ = The moments m s satisfy the infinite system of linear 

equations 

m s = ) \{-\) L c L+s \~ ' " " \m L , s>l. 



-a-tf-f::: 1 ) 
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Proof. Indeed, for %z < we have (recall that m = 1): 



- G(z) = Z ™ — t + X ™ — — ) = Z ^ Z m 4 — ) L+1 - 

n=\ n=l «=1 L=0 

This series is absolutely and uniformly convergent for !Rz < 0, as is implied by (PTOl) . We obtain the 
needed result after taking the sth left derivative at z = 0. ■ 

Numerical calculations are presented in J3]|. This method gives high precision values for other 
constants, including the Kinney's constant. 

5. Exponential generating function m(t) 

The aim of this section is to interpret ([141) in terms of m(t). The following Theorem, along with 
the boundary condition m(0) = 1 and regularity property as in Lemma [Q, uniquely determines the 
function m(t) . 

Theorem 2. The function 111(5) satisfies the integral equation 

m(-s) = {2e s -1) J m' (-0/o(2 yfsl)dt, s e R+, (16) 


where Jq(*) stands for the Bessel function: Jq(z) = \ cos(z sin x) dx. 
Proo/ For %z < 1, we have that G(z) = J^°° m\zi)e~ 1 dt. Thus, 

00 00 

G(z) = -- f m'(-?y /z di for Slz < 0, G(z) = - f m'(f)e~' /z dffor < %z < 1. 
z J z J 



Thus, the functional equation for G(z) in the region %z < -1 in terms of m'{t) reads as 

00 

1 C 2 11 

- = m'(-0f -e^ + -e tz - -e'Adt. (17) 

z J v z + 1 z z 7 



Now, multiply this by e"* z and integrate over %z = —cr < -1, where 5 > is real. We have ([30|, p. 
465) 



—cr-tioo 

19 



dz = -2ni\ 



cr-l+ioo <T +ico 



2 J -— pdz = -2e' s I — dz = -2e s I — dz = -4me s J (2 Vtf), 



0--l-;oo (TQ-ioo 
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where cr = {a - I) > 0, and Ja(*) stands for the Bessel function (see POL p. 597 for the 
representation of the Bessel function by this integral). Further, 



-sz 



■dz = -2niJ (2yfsi). 



-cr+ioo ( -cr+ioo 

r e (t-s) z ^ I _27n if s > t, C 
J z |o if s<t, J ■ 

—<T—ioo V —cr—ioo 

Thus, eventually 

oo s 

-2m = -2m J m'(-t)(2e s - 1)7 (2 V7t) dt - 2ni J m'(-t) dt; 
o o 
since m(0) = 1, this proves the Proposition. ■ 

Thus, we have obtained an integral equation for m(s), which corresponds to the functional equation 
(fl4l) for G(z). They are both, in fact, mere reformulations of (Q~|) in different terms. 

6. Dyadic eigenfunctions 

In this Section we introduce the sequence of functions G A (z), which satisfy the functional equation 
analogous to (PT41) . 

Since J' Q (*) = -Ji(*), Ji(0) = 0, integration by parts in (PT6l) leads to 

CO 



J vt vs 



(18) 



yfs Vs(2e s - 1) 
o 

J "•CO /— 
f(r)yrdr 

converges absolutely) is defined as 

oo 

gip) = J f(r)J v (rp)rdr, 
o 

where 7 V (*) stands for the vth Bessel function. The inverse is given by the Hankel inversion for- 
mula with exactly the same kernel ([49], chapter XIV, section 14.4.). Thus, after a proper change of 
variables, Hankel transform reads as 

CO CO 

g(p) = J f(r)J v (2yfw)dr^f(r) = J g(p)J v (2 ^) dp. 
o o 
Thus, application of this inversion to the identity (fT8~l) yields 

oo oo 

V* J a/7 J V?(2e'-1) 

o o 
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The first integral on the r.h.s. is equal t —^J Q (2y[si) |~ = 4=. Let t//(s) = (2e s - 1) 1/2 . Then this 
equation can be rewritten as 

oo 

m(-s) 1 f m(-f) 7i(2VH) , 



f m(- 

J 



Hence, if we denote 



o 



7i(2V^) % m(-j) - 1 

£O,0, — - - - = m(s), 



iKs)*ff(t) V#0) 
we obtain a second type Fredholm integral equation with symmetric kernel ( ||2"5ll , chapter 9): 



CO 

-I 



m(s) = £(s) - m(t)K(s, t) dt, 



o 



where £(s) = 



oo 

i r *eV5) , ■ ( y,- lr .,). 



The behavior at infinity of the Bessel function is given by an asymptotic formula 

/ 2 xl/2 / 3 v 

Ji(x)~[— ) cos(x--Tr) 
(EH, chapter VII, section 7.1). Therefore, obviously, 

OO OO oo 

J J \K(s,t)\ 2 dsdt < oo, J~ |^)| 2 ds<oo. 



Thus, the operator associated with the kernel K(s,t) is the Hilbert-Schmidt operator ( [|25l , p. 532). 
The theorem of Hilbert-Schmidt ( Il25ll . p. 283) states that the solution of this type of integral equations 
reduces to finding the eigenvalues A and the eigenfunctions A A (s). We postpone the solution of this 
integral equation for the future. Till the end of this section we deal only with eigenfunctions. The 
integral operator, consequently, is a compact self-conjugate operator in the Hilbert space, it possesses 
a complete orthogonal system of eigenfunctions A, all A are real and A n — » 0, as n — » oo. If we denote 
A A {s)if/{s) = B A (s), then the equation for an eigenfunction reads as 



CO 

I 



2e' - 1 



This gives 5^(0) = 0. Since A A (s) e L 2 (0, oo), and 7i(*) is bounded, this implies that B A (s) is 
uniformly bounded for s > as well. Moreover, since the Taylor expansion of 7i(*) contains only 
odd powers of the variable, B A (s) sfs has a Taylor expansion with center and is an entire function-h 
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Now, multiply this by yfse s/z , z > 0, and integrate over s e R + . The Laplace transform of yfs J i(2 yfs) 
is j?e~ 1/z (|(30l, p. 503). Thus, we obtain 

CO CO 

l/t^-' ; *4j^>V^d, (19) 



Denote by G A (-z) the function on both sides of this equality. Thus, G A {z) is defined at least for 
%z < 0. Since 2e t{z+l) - e tz = (2e' - l)e tz , we have 

CO 

A(2G A (z + 1) - G A (z)) = J B A (t) <te tz dt = ^G x (l/z). 

o 

Therefore, we have proved the first part of the following Theorem. 

Theorem 3. For every eigenvalue A of the integral operator, associated with the kernel K(s, t), there 
exists at least one holomorphic function G A (z) (defined for z e C \ R>J, such that the following holds: 

2G A (z + 1) = G A (z) + i<?4-)- (20) 

Az l x z' 

Moreover, G A (z)for %z < satisfies all regularity conditions imposed by it being an image under the 
Laplace transform ( If30~1 , p. 469). 

Conversely: for every A, such that there exists a non-zero function, which satisfies (|20l) and these 
conditions, A is the eigenvalue of this operator. The set of all possible A's is countable, and A n — > 0, 
as n — > oo. 

Proof. The converse is straightforward, since, by the requirement, G A (z) for %z < is a Laplace 
image of a certain function, and all the above transformations are invertible. We leave the details. If 
the eigenvalue has multiplicity higher than 1, then these /l-forms span a finite dimensional C-vector 
space. Note that the proof of Proposition Q] implies \A\ < 0.342014 + . Finally, the functional equation 
(|20l) gives the analytic continuation of G A (z) to the half -plane %z < 1. Further, if z 6 1A, where 
tl = {0 < %z < 1} \ {|z| < 1}, we can continue G A (z) to the region 14 + 1, and, inductively, to 14 + n, 
n 6 N. Let 14 be the union of these. We can, obviously, continue G A (z) to the set 14 Q 1 + n, n e N. 
Similar iterations cover the described domain. ■ 

Note that, in contrast to G(z), we do not have a symmetry property for G A {z). 

Next calculations produce the first few eigenvalues. Let the Taylor expansion of G A (z) is given by 

oo 
L=l 

It converges in the unit circle, including its boundary (as is clear from (T20l) . there exist all left deriva- 
tives at z = 1). Thus, m^ have the same vanishing properties as m L (which guarantees the convergence 
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of the series in (flOl)). And therefore, as in the Proposition [51 we obtain: 

Amf^p-l^c^y/^y?, s>l. 

Here c L = 2,T=i zt?- ^ we denote e StL = (-l)^ 1 ^^**" 1 ), then A is the eigenvalue of the infinite 
matrix &™ L=l - The numerical calculations with the augmentation of this matrix at sufficiently high 
level give the following first eigenvalues in decreasing order, with all digits exact: 

A 1 = 0.25553210+, A 2 = -0.08892666+, A 3 = 0.03261586+, A 4 = -0.01217621+. 

Paper [1 J contains graphs of G A (z) for the first six eigenvalues, and, more importantly, "pair-correlation" 
results among different eigenvalues, including the "eigenvalue" -1, which corresponds precisely to 
G{z). These results reveal the importance of Ga(z) in the study of F(x). 

7. p-ADIC DISTRIBUTION 

In the previous sections, we were interested in the distribution of the nth generation of the tree T 
in the field of real numbers. Since the set of non-equivalent valuations of Q contains a valuation asso- 
ciated with any prime number p, it is natural to consider the distribution of the set of each generation 
in the field of p-adic numbers Q p . In this case we have an ultrametric inequality, which implies that 
two circles are either co-centric or do not intersect. We define 

F n (z, v) = 2"" +1 #{ a - 6 T (n) : ord p (^ - z) > v}, z e Q P , v e Z. 
b b 

(When p is fixed, the subscript p in F n is omitted). Note that in order to calculate F n (z, v) we can 
confine to the case ord p (z) < v; otherwise ord p (| - z) > v <=> ord p (|) > v. We shall calculate the limit 
distribution p. p {z, v) = lim,,^ F n (z, v), and also some characteristics of it, e.g. the zeta function 

Z p (s) = J \u\ s d/u p , seC, ze Q p , 

ueQ p 

where | * | stands for the p-adic valuation. 

To illustrate how the method works, we will calculate the value of F n in two special cases. Let p = 2 
and let E(n) be the number of rational numbers in the nth generation with one of a or b being even, 
and let 0(n) be the corresponding number fractions with both a and b odd. Then E(n) + 0(n) = 2"" 1 . 
Since | in the nth generation generates ^ and ^ in the (n + l)st generation, each fraction | with 
one of the a, b even will generate one fraction with both numerator and denominator odd. If both a, 
b are odd, then their two offsprings will not be of this kind. Therefore, 0(n + 1) = E{n). Similarly, 
E(n + 1) = E(n) + 20(n). This gives the recurrence E(n + 1) = E(n) + 2E(n - 1), n > 2, and this 
implies 

2" + 2(-l)" 2"' 1 + 2(-l)"- 1 2 
E(n) = 0(n) = 3 V // 2 (0,0) = -. 
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(For the last equality note that | and - simultaneously belong to T (n \ and so the number of fractions 
with ord 2 (*) > is E(n)/2). We will generalize this example to odd prime p > 3. Let L,(n) be the 
part of the fractions in the nth generations such that ab~ l = i mod p for < i < p — 1 or i = oo (that 
is, b = mod p). Thus, 

ieFpUoo 

in other words, L,(n) = F n (i, 1). For our later investigations we need a result from the theory of finite 
Markov chains. 

Lemma 3. Let A be a matrix of a finite Markov chain with s stages. That is, > 0, and Yfj=\ a i,j - 1 
for all i. Suppose that A is irreducible (for all pairs (i, j), and some m, the entry ajj of the matrix 
A'" is strictly positive), acyclic and recurrent (this is satisfied, if all entries ofA m are strictly positive 
for some m). Then the eigenvalue 1 is simple and if A is another eigenvalue, then \A\ < 1, and A" 1 , 
as m — * oo, tends to the matrix B, with entries bij = 7ij, where (tti, ...,n s ) is a unique left eigenvector 
with eigenvalue 1, such that Yfi=\ n j = 1- 

A proof of this lemma can be found in ll20ll . Section 3.1., Theorem 1.3. 

Proposition 6. p p (z, 1) = -p^for z 6 Z p . 

Proof. Similarly as in the above example, a fraction | from the nth generation generates -^r and 
%Y in the (n + l)st generation, and it is routine to check that 

L i (n + l) = ^L i t,(n) + ^L i _ l (n) for i e F, U {«,}, (21) 

(Here we make a natural convention for j^. and i - 1, if i = 1 or oo). In this equation, it can happen 
that i - 1 = y~: mod p; thus, (2i - l) 2 = -3 mod p. The recurrence for this particular i is to be 
understood in the obvious way, L,(n + 1) = L,_i(n). Therefore, if we denote the vector-column 
(Loo(n), L (n), L p -i(n)) T by v„, and if Jl is a matrix of the system (f2T|) . then v„+i = Jl\ n , and hence 

where Vi = (0, 0, 1, 0, 0) r . In any particular case, this allows us two find the values of L, explicitly. 
For example, if p = 7, the characteristic polynomial is 

f(x) = -L(jc - l)(2x - \)(2x 2 + \)(4x 4 + 2x 3 + 2x+ 1). 
16 

The list of roots is 

1 i -1 - VF7 V 1 + ^ 

2 V2 8 2V2 
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(with respect to the two values for the root Vl7), the matrix is diagonalisible, and the Jordan normal 
form gives the expression 

8 

Uin) = V Cuati. 



=1 



Note that the elements in each row of the (p + 1) x (p + 1) matrix Jl are non-negative and sum up 
to 1, and thus, we have a matrix of a finite Markov chain. We need to check that it is acyclic. Let 
t<7) = i - 1, and = for i e ¥ p U {oo}. The entry aff ofW is 

(m) _ V"< 

'\, — ,'m-\ 

Therefore, we need to check that for some fixed m, the composition of m cr's or t's leads from any i 
to any j. One checks directly that for any positive k, and i, j G ¥ p , 

/"^oo-OT^ocroT''" 1 © = j, 
T p ~~ 1 ~ j ocro /(oo) = j, 
r octot'" 1 ® = oo; 

(for i = 0, we write r" 1 for r''" 1 ). For each pair (i, j), choose k in order the amount of compositions 
used to be equal (say, to m). Then obviously all entries of J?T" are positive, ant this matrix satisfies the 
conditions of Lemma[3l Since all columns also sum up to 1, (tti, ...,n p+ i), nj = j-^, 1 < j < p + 1, is 
the needed eigenvector. This proves the Proposition. ■ 

Next Theorem describes fi(z, v) in all cases. 

Theorem 4. Let veZ and z e Q p , and ord p (z) < v (or z = 0). Then, ifz is p-adic integer, 

1 

Kz, v) = 



p v + pv-l 



Ifz is not integer, ord p {z) = -A <0, 

Viz, v) = 

For z = 0, -v < 0, we have 



pV+2A _|_ pV+2A-\ 



M0,-v) = i- * - . 

pv+l + pV 



This Theorem allows the computation of the associated zeta-function: 
Corollary 1. For s in the strip -1 < %s < 1, 

andZ„(s) = Z p (-s). 



22 



GIEDRIUS ALKAUSKAS 



The proof is straightforward. It should be noted that this expression encodes all the values of p(0, v) 
for veZ. 

Proof of Theorem^ For shortness, when p is fixed, denote ord p (*) by v(*). As before, we want 
a recurrence relation among the numbers F n (i, k), i e Q+. For each integral k, we can confine to the 
case i < p K . If i = 0, we only consider k > and call these pairs (i, k) "admissible". We also include 
G„(0, -k) for k > 1, where these values are defined in the same manner as F n , only inverting the 
inequality, considering | e T ( "\ such that v(|) < —k; this is the ratio of fractions in the nth generation 
outside this circle. As before, a fraction | in the nth generation generates the fractions and ^ in 
the (n + l)st generation. Let r(z, k) = ((/ - 1) mod p K , k). Then for all admissible pairs (i, k), i ± 0, the 
pair r(i, k) is also admissible, and 

v(— i) = /c v(- - (i - 1)) = /c. 

Second, if = i + p K u, i+ I, u e Z p , and (7, /c) is admissible, then 

a i p"u 

b p K u) 

Since v( ^) = v(z) - v(l - z), this is unless z is an integer, equals to v(z) if the latter is > and equals 
to -v(l - z) if v(l - z) > 0. Further, this difference has valuation > kq = k, if i e Z,i S 1 mod p, 
valuation > kq = k - 2v(l - z), if i e Z, i = 1 mod p, and valuation > a: = k - 2v(z) if i is not 
integer. In all three cases, easy to check, that, if we define z = -j-r- m od /j* , the pair <x(z, a:) = def (z'o, k ) 
is admissible. For the converse, let | = z + /j*°m, m 6 Z p . Then 

a z' p K ° 



a + b 1+z'o (1 + z'o + p Ko u){\ + z' ) 
If i = is a p-adic integer, z ^ 1 mod p, this has a valuation > a: = /c ; if z is a p-adic integer, 
z = l(p), this has valuation > k = k - 2v(z ) = ato + 2v(l - z); if z is not a p-adic integer, this has 
valuation > k = k - 2v(l + z ) = Kq + 2v(z). Thus, 

V( - l) > K V(- - Z ) > AT . 

a + b b 

Let z = 1. If = 1 + p K zi, then /c > 0, u 6 Z,„ and we obtain | = -1 — < -/<-. Converse is 
also true. Finally, for /c > 1, 



and 



a + b a 
v(— — ) < -/c « v(-) < -/c, 
o b 



a a 

v( -) < -k & v(- + 1) > AT. 

a + /> o 



Therefore, we have the recurrence relations: 

' F n+1 (i, k) = jF n (r(i, k)) + jF n (cr(i, k)), if (z, *-) is admissible, 
■ F n+1 (l,K) =\F n {0,K) + \G n {0,-K\K>\, (22) 

k G n+ y{Q,-K)=\G n {Q,-K) + \F n {-\,K),K>\. 
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Thus, we have an infinite matrix 3\, which is a change matrix for the Markov chain. If v„ is an infinite 
vector-column of F' n s and G' n s, then v n+ i = Jl\ n , and, as before, v„ = 3i n ~ x \\. It is direct to check that 
each column also contains exactly two nonzero entries |, or one entry, equal to 1. In terms of Markov 
chains, we need to determine the classes of orbits. Then in proper rearranging, the matrix J{ looks 
like 



(Pi 


. 


. o ..: 





p 2 


. ... 





. 


• P.v 


v • 




• 



where P 4 are finite Markov matrices. Thus, we claim that the length of each orbit is finite, every orbit 
has a representative G*(0, -k), k > 1, the length of it is p K + p K ~ l , and the matrix is recurrent (that 
is, every two positions communicate). In fact, from the system above and form the expression of the 
maps r(7, k) and cr(i, k), the direct check shows that the complete list of the orbit of G*(0, -k) consists 
of (and each pair of states are communicating): 

G*(0, -k), 
F*(i,K) (i = 0,l,2,...,p*-l), 
F*(p~ a u,k-2A) (A = 1,2,..., a:- l,u e N,u £ mod p,u < p K ' A ). 

In total, we have 

K-i 

i+ p k + Y}-p k ~ a - p k ~ a ~ 1 ) =p k+ p k ~ 1 

members in the orbit. Thus, each P„ in the matrix above is a finite dimensional i K x l K matrix, where 
4 = p* + p K ~ x ■ For k = 1, the matrix P\ is exactly the matrix of the system (I2TI) . As noted above, the 
vector column (1, 1, l) T is the left eigenvector. As in the previous Proposition, it is straightforward 
to check that this matrix is irreducible and acyclic (that is, the entries of P" are strictly positive for 
sufficiently large n). In fact, since by our observation, each two members in the orbit communicate, 
and since we have a move G t (0, -k) — » G*(0, -k), the proof of the last statement is immediate: there 
exists n such that any position is reachable from another in exactly n moves, and this can be achieved 
at the expense of the move just described. Therefore, all entries of P" are strictly positive. Thus, the 
claim of the Theorem follows from the Lemma |3j ■ 

8. Conclusion 

We end the paper with the following remarks. As is implied by TheoremUl the measure fi p of those 
rationals in the Calkin-Wilf tree which are invertible elements of Z p is equal to We follow the 
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line of the Tate thesis J9]|, and modify this measure in order Z* to have measure 1; accordingly, let us 
define fi' p = jr x li P - Thus, we are lead to the formal definition of the zeta function 

Ms) - n L r ■ n (• - ?) n ■ iz^r = > + w- + 

This product diverges everywhere; nevertheless, if we apply the functional equation of the Riemann 
f function for the second multiplier, we obtain 

frW = ^(2^r $ cos(y)r(^(^(5 + 1). 

From the above definition it is clear that, formally, this zeta function is the sum of the form 2reQ+ /V~* \ 
where, if r e Q + , and fi r stands for the limit measure of those rationals in the nth generation of T, 
which have precisely the same valuation as r at every prime which appears in the decomposition of 
r, times the factor riord p (r)#o ~[- Surprisingly, the product £(s)g(s + 1) is the zeta function of the 
Eisenstein series G\{z), which is related to the distribution of rationals in T at the infinite prime 
Qeo = R. In fact, 

CO 

J (Gi(fe) - G^ioo)^- 1 dz = -87r 2 (2;r)-T<>)<r<>)^> + 1). 
o 

This is a strong motivation to investigate the tree T and the Minkowski question mark function in a 
more general - idelic - setting, thus revealing the true connection between p-adic and real distribu- 
tion, and clarifying the nature of continued fractions in this direction. We hope to implement this in 
the subsequent papers. 

Unfortunately, currently we left the most interesting question (the explicit description of the mo- 
ments of ?(jc)) unanswered. It is desirable to give the function G(z) and dyadic forms Ga(z) certain 
other description than the one which arises directly from the tree T . This is in part accomplished 
in [0Q| and 0. These papers form a direct continuation of the current work. Among the other re- 
suits, the dyadic zeta function 1S introduced: it is given by £m( s W(s + 1) = J x s dF(x); the 

nature of dyadic eigenfunctions Gx(z) is clarified; certain integrals which involve F(x) are computed; 
and finally, this research culminates with the proof that in the half plane %z < 1 the dyadic period 
function G(z) can be represented as an absolutely convergent series of rational functions with rational 
coefficients. Possibly, this technique can find its applications in the study of period function for Maass 
wave forms. 
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